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Abstract 

Anisotropic lattices, with a temporal lattice spacing smaller than the spatial one, 
allow precision Monte Carlo calculations of problems that are difficult to study 
otherwise: heavy quarks, glueballs, hybrids, and high temperature thermodynamics, 
for example. We here perform the first step required for such studies with the 
(quenched) Wilson gauge action, namely, the determination of the renormalized 
anisotropy £ as a function of the bare anisotropy £o and the coupling. By, essentially, 
comparing the finite-volume heavy quark potential where the quarks are separated 
along a spatial direction with that where they are separated along the time direction, 
we determine the relation between £ and £o to a fraction of 1% for weak and to 
1% for strong coupling. We present a simple parameterization of this relation for 
1 < £ < 6 and 5.5 < (3 < oo, which incorporates the known one-loop result and 
reproduces our non-perturbative determinations within errors. Besides solving the 
problem of how to choose the bare anisotropies if one wants to take the continuum 
limit at fixed renormalized anisotropy, this parameterization also yields accurate 
estimates of the derivative <9£o/<?£ needed in thermodynamic studies. 



1 Introduction 

Lattice QCD has experienced rapid progress in the last few years (see e.g. the proceedings 
of the latest lattice conference Q). In particular, it should be viewed as progress that we 
can now clearly see the limitations of various approximations that have been commonly 
used in the past, like quenching or the use of (mean-field improved) perturbative 
coefficients for the improvement terms in actions and currents. 

A great theoretical advance in the last two years was the determination of the non- 
perturbatively 0(a) improved Wilson quark action for both quenched (2], |3|, Q and full 
(nf = 2) QCD H initiated by the ALPHA collaboration. These results are most relevant 
for the study of light hadrons. They are therefore only of limited value in the study 
of heavy quark systems, where we are also beginning to see the limitations of current 
methods. The most popular of these are the NRQCD Q and Fermilab @ approaches. 



The problem is, as for the standard Wilson quark action, that the non-perturbative 
corrections to the 0(a) improvement coefficients in these actions are large on coarse 
lattices. Different mean-field prescriptions to take some of these corrections into account 
give quite different results on coarse lattices §. A non-perturbative determination of 
these coefficients is therefore called for. 

We think that the most reliable and ultimately simplest approach is to perform the 
non-perturbative 0(a) improvement not for the NRQCD or Fermilab formalisms, but 
rather to extend the Symanzik improvement program for Wilson-type quark actions to 
anisotropic lattices, with a temporal lattice spacing at much smaller than the spatial one 
a s . On such lattices heavy quarks will not suffer large lattice artifacts as long as their 
masses are small in units of a^. At the expense of the relatively modest cost incured by 
simulating a lattice with a larger temporal extent one can study heavy quark systems in a 
relativistic framework with controlled errors. Such studies would be orders of magnitude 
more expensive on isotropic lattices. The use of anisotropic lattices for heavy quark 



physics was advocated in refs. [g, 1C, 11]; we refer to || for studies of the classical case 
and to (T(| for some exploratory simulations. 

Anisotropic lattices are also important in other situations. Generally speaking, they 
allow one to reap many of the benefits of fine lattices, while still using cheap coarse 
spatial lattices: 

• The signal to noise ratio of correlation functions calculated in Monte Carlo 
simulations generically decays exponentially in time. Choosing a smaller temporal 
lattice spacing gives more time slices with an accurate signal, enabling more precise 
and confident mass determinations, for example. This is crucial for particles with 



bad signal/noise properties like glueballs [12], hybrids and P-state mesons. For the 



same reason it is also clear that excited-state masses can be determined much more 
accurately on anisotropic lattices. 

• For a full determination of thermodynamic quantities one has to be able to take 
independent derivatives with respect to temperature and volume. The simplest 
way to achieve this, is to have independent spatial and temporal lattice spacings. 
The point of the previous paragraph also becomes important at high temperatures, 
where one needs a small temporal lattice spacing to have a sufficient number of 
Matsubara frequencies for accurate studies of the electro-weak phase transition, 
transport coefficients in the quark-gluon plasma phase, and many other issues. 

• By euclidean invariance we can also think of the small lattice spacing as a spatial 
one. This allows one to study large spatial momenta, which are phenomenologically 
important for form factors, for example. 

Of course, the advantages of anisotropic lattices come at a price, namely, that 
coefficients in the action have to be tuned to restore space-time exchange symmetry on 
the quantum level. In the context of the 0(a) improvement program, the good news is 
that the increase in the number of coefficients that have to be determined for a relativistic 
action is quite modest: For a Wilson- type quark action there are three (instead of one) 
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coefficients to be tuned to achieve non-perturbative 0(a) improvement [11|. How to 



achieve this using the Schrodinger functional, similar to the isotropic case H, has been 



outlined in pl|. 



However, before worrying about the quark action, we have to consider gauge actions 
on anisotropic lattices. The first step is the study of pure gauge theory, relevant for 
quenched QCD. In this paper we will consider the simplest case, the anisotropic Wilson 
gauge action [14|].[] For this action no coefficients have to be tuned to restore space- 
time exchange symmetry up to 0(a 2 ) errors. Nevertheless, there is something to be 
done, since we have to know the true or renormalized anisotropy £ = a s /at as a function 
of the bare parameters. This is important, because we would like to take the continuum 
limit at a fixed value of the renormalized anisotropy, for example. 

Our aim in this paper is to develop a simple and accurate method to determine the 
renormalized anisotropy and apply it to the Wilson plaquette action. For SU(iV) gauge 
fields we write this action on an anisotropic lattice as 

S = i E ^ReTr(l-P ss ,(x)) + £oReTr(l-P 0s OE)). (1.1) 

N x,s>s> W 

Here P^(x) is the plaquette operator in the ^-plane.Q i.e. the product of link fields 
around a single plaquette (with x chosen according to any consistent convention as one 
of its corners). Expanding the action in terms of the field strength for small lattice 
spacings, it is easy to see that at the classical level £o = &s/ a *j ^ we use P = 2-/V/g 2 
to introduce the bare coupling g as in the isotropic case. This justifies the name bare 
anisotropy for £o- 

Our problem is to determine the relation £ = £(£o,/3), or its inverse, £o = £q(£,(3), 
which is actually preferrable from a practical point of view. In sect. 2 we describe 
the physical ideas behind the "ratio method" we use to determine this relation non- 
perturbatively. Basically the same method was used previously in [l^] (but not 

always with the same results, as we will see). In sect. 3 we present the details of our 
implementation of these ideas and of our simulations of the anisotropic SU(3) plaquette 
action. Our results cover both fine and coarse lattices and five anisotropies in the range 
1-5 < £ < 6. In sect. 4 we use known perturbative results to motivate paramaterizations 
that reproduce all our data within errors. We also briefly discuss the application of our 
results in high precision thermodynamic studies. We conclude in sect. 5 with a summary 
and an outlook on future work. 



1 Determinations of the renormalized anisotropy for improved actions will be reported in |l5| . 

2 Here a summary of our lattice conventions: We always work on four-dimensional hypercubic lattices 
with periodic boundary conditions. Directions are labelled by ft — 0, 1, 2, 3 or t, s, s' , where t = stands 
for the time and s, s' for a space direction. a M is the lattice spacing in direction fx (we write if we want 
to consider the spacing as a vector in the positive jj, direction), and we assume spatial isotropy, i.e. that 
all spatial a M are the same, equal to a s . For the extensions of the lattice we write L M and iV M = i M /a M 
in physical and lattice units, respectively. 
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2 Non- Pert ur bat ive Determination of the Renormalized 
Anisotropy 



The obvious spectral quantity to extract the renormalized anisotropy £ of a gauge action 
from is the static potential. Actually, on an anisotropic lattice there are two potentials, 
according to whether the heavy quark and anti-quark propagate along the time or a 
space direction. We will refer to them as the "regular" and "sideways" potential, and 
denote them by Vt (r) and V s (r), respectively. 

On the lattice Vf(r) is measured in units of the temporal lattice spacing and V^(r) 
in terms of the spatial one. They therefore differ by a factor of £. However, they also 
differ by an additive constant, since the self-energy corrections to the static potential 
are different if the quarks propagate along the time or a space direction. To determine £ 
from a comparison of the regular and sideways potential we therefore have to demand, 

V^r) = £Vt(r) + const , for sufficiently large spatial r . (2-1) 

Note that due to the 0(a^,af) errors of the potential(s), the above equality will in 
principle only hold for asymptotically large r. From experience with the static potential 
we expect that the components of r have to be at least 2 — 3 in units of a s , for the 
systematic 0(a%, af) errors in £ to be below the 1% level. 

Due to the constant in fl2.1|) this approach always requires several r values to obtain 
an estimate of £. Since one cannot use small r and errors increase rapidly at larger r, this 
approach is relatively expensive (and not very elegant, with a complicated error analysis 
if done properly). 

A much better approach involves comparing the sideways potential with the quarks 
separated along a spatial direction with the case where they are separated along the time 
direction.^ There is no additive constant to complicate matters and we can demand 

V s (y&s) = Vs{y£,&t) for sufficiently large y . (2.2) 

Using one's favorite method of calculating the static potential one can obtain an estimate 
of £ for each value of y (on the right hand side one will generically have to interpolate 
between two values of y£). The estimated £ should rapidly reach a plateau as y increases. 

This method is perfectly fine, in principle and in practice. If V, is strictly meant to 
be the (sideways) static potential with negligible finite-volume effects, it does however 
have a limitation. Namely, it will be difficult to reach small lattice spacingsj^j where the 
elimination of finite-volume effects in the potential becomes very expensive. However, 
the finite-volume potential is also a spectral quantity; if we can arrange the physical 
extent of the lattice to be the same in the spatial and the time directions, then the 
finite- volume effects in V s (ya s ) and V s (y^at) should also be the same! (Below we will 

3 We could also allow the separation to have a component along (another) spatial direction. 

4 As briefly mentioned in the introduction, we would like to reach small lattice spacings and make 
contact with perturbation theory, so that we can present an analytic formula for the relation between £ 
and £o from weak to strong coupling. 
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describe one way of getting around the problem of how to choose the size of the lattice to 
correspond to the same physical extent in all directions without knowing £ beforehand.) 

So, finite- volume effects should not prevent as from measuring £ on fine lattices. The 
remaining question is if we have a signal for "times" large enough so that ratios of Wilson 
loops allow us to extract the asymptotic value of the (finite-volume) potential. Actually, 
this is not necessary. Similarly to the logic that allowed us to work with the finite- volume 
potential we note the following. If the ratios of Wilson loops that asymptotically would 
reach V s (ya s ) and V s {y^at) have the same excited-state (i.e. finite- "time" ) contributions, 
we can use these ratios at finite "times" to measure £. In other words, we have to set up 
a situation where the heavy quark anti-quark states that can be thought of as underlying 
the spatial and temporal Wilson loops in question correspond to the same physical state 
from a continuum point of view, i.e. when ignoring the usual 0(a1,af) lattice artifacts. 

Before proceeding, let us introduce some notation for ratios of Wilson loops, namely, 

W ss (x+l,y) W st (x+l,t) 

in terms of spatial, W ss , and temporal, W s t, Wilson loops. To avoid cumbersome 
notation we here measure x,y in units of a s and t in units of at- Asymptotically, for 
large x, the ratios R ss (x, y) and R s t(x, t) approach exp[—a s V s (ya s )] and exp[— a s V s (t&t)], 
respectively. 

For R ss and R s t to have the same excited-state contribution the spatial and temporal 
links making up the Wilson loops in ( [2 .3D have to be smeared by the same amount 
in physical units. A simple way of achieving this is with no smearing at all. This 
leads to the "ratio-of- Wilson-loop" method for determining £, which has been used 
previously |16, [P7|, |18]. 



We will describe our detailed implementation(s) of this method in the next section, 
where we also provide empirial support for our statements concerning the cancellation 
of finite "time" and volume effects in R ss (x,y) and R s t(x,t). 



3 Simulation and Results 

Let us start by summarizing our discussion of the ratio method in the previous section. 
Although the ratios of Wilson loops, R ss (x,y) and R s t(x,t) of eq. (fO|), are not spectral 
quantities (for x where they can be measured reasonably accurately in practice), we 
expect excited-state corrections to cancel between them if t = £y and all gauge links are 
unsmeared. Similarly, finite-volume corrections to R ss (x,y) and R s t(x,t) are the same 
if the temporal and spatial extents of the box are equal in physical units, i.e. Nt = £iV s 
in lattice units. Of course, these statements are expected to hold only if x, y and t are 
not too small; otherwise there can be large 0(a^,a^) lattice artifacts. 

The conceptually cleanest way to implement the above ideas in a simulation is to 
turn the problem around and not measure £ for fixed £o> but rather try to determine the 
£o that corresponds to a given £. The following method is simplest if £ is an integer, so 
the reader might want to have this case in mind for the moment. 
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• For fixed (3 and £ choose a lattice volume N% x iV* with N t = £,N s and calculate 

5(x,y) = <5(x,y|£) = 8(x,y;t = £y) , where S(x,y;t) = ^^'^ - 1, 

R st {x,t) 

(3.1) 

for two or three trial values of £o- 

• Interpolate in £o to find the zero-crossing of S(x,y) for fixed x,y. This determines 
an estimate of the non-perturbative £o(£>/?) for the given 

• Consider different x,y; look for a plateau of £o(£>/3) as x and y increase. 

By starting at weak coupling, where perturbation theory can be used, and working 
one's way towards large coupling it is relatively easy to estimate trial values of £o that 
bracket the correct £o(£, f3) quite closely. In almost all cases we only needed two trial 
values for £o- The cases where simulations were performed at more than two values were 
used to check that the 5{x,y) are sufficiently linear as functions of £o over the regions 
we typically had to interpolate. 

We are mainly interested in integer £ for which the above method is best suited. 
However, there is no problem in using it also for "small fractions". For example, one 
case we investigated somewhat is £ = §. For even N s almost everything goes through 
as above; the only change is that for S(x,y = 3), for example, one has to perform an 
interpolation of 5(x, y = 3; t = 4) and 5(x, y = 3; t = 5) to i = £-y = 4.5. 

As remarked earlier, the above method is perhaps the conceptually cleanest way to 
proceed, since we choose £ not £o from the start and can therefore make sure that L s = L t 
in physical units. It was our method of choice for all but the coarsest lattices. However, 
on the coarsest lattices, where one needs at least 1000 configurations to measure £o(£) 
with some accuracy, one would like to avoid having to do independent simulations at two 
or more trial values of £o (on fine lattices one only needs a few hundred configurations, 
so this is not such an important issue). This can be achieved with a slightly modified 
procedure: 

• Given a good guess for the value of £o corresponding to the desired £, calculate 
R ss (x,y) and R s t(x,t) by simulation at £o on an iV s 3 x N t lattice with N t = £N s . 

• For given x,y interpolate R s t(x,t) in t to match R ss (x,y). Estimate £ as £ = t/y. 

• Find the plateau of the estimated £ as x and y increase. 

As long as our estimate of £o was n °t too bad, all our previous remarks about the 
cancellation of finite volume effects and such should apply. In this method we obtain 
£(£o 5 0), which can be translated into an estimate of £o(£j P) at the desired nearby £. We 
have explicitly checked in a few cases that both methods give consistent results. 

We performed simulations at a large set of couplings aiming at £ = 2,3,4 and a 
smaller set for £ = 1.5, 6. We generated between a few hundred and up to 2700 (almost) 
independent configurations at weak, respectively, strong coupling. We used two sets of 



6 



1.05 



1 1 


1 1 1 1 1 1 1 1 1 1 
18 = 18.0, 16 3 x32, f = 1.939 


1.10 




i i i i i i i i i i 

(3=18.0, 16 3 x32, f =1.939 




x R ss (x,3)/1.13 




x R ss (x,4)/1.13 










— st, \ > J / * 


- g 


o R gs (x,3)/R st (x,6) 

H 

H h m H H 


1.05 


- 


o R ss (x,4)/R st (x,8) - 

N 

H H H H _ 





- -o e & e -0 — 

i i i i i i i i i i 


1.00 


o 


- - -Q & B -0 

ill l l l l l l l l l 



2 4 6 2 4 6 

X X 



Figure 1: Ratios R ss (x,y), R s t(x,t) and the "ratios of ratios" R ss (x,y)/ 'R s t(x,t = £y) as 
functions of x for fixed y,t for a very fine lattice, f3 = 18, £ ~ 2. Shown are y = 3 (left) 
and y = 4 (right). The R ss (x,y) and R s t(x,t) have been rescaled by some overall factor 
simply so that they can be put more meaningfully on the same plot with their ratio. 



code. One employs Metropolis, the other Kennedy- Pendleton IS] heatbath updating; in 
both cases alternating with micro canonical over-relaxation steps 2C]. 

We started at weak coupling where our method worked so well that we initially did 
not use link integration pl| for Wilson loops. On coarse lattices, however, fluctuations 
are large and link integration yields a large saving in CPU time. It was implemented by 
replacing independent links in a Wilson loop by their average over 20 local updates at 
fixed values of the staples surrounding the given link. 

Most simulations were performed on lattices with N s = 8 or 10. In several cases we 
also ran at N s = 12, 16, with otherwise the same parameters, and found no significant 
difference from the results on smaller volumes. This is in accord with the theoretical 
expectations discussed earlier. 

We performed one fine, one intermediate, and one coarse lattice simulation with 
N s = 16, generating about the same number of configurations as for the smaller volumes 
(we used three trial values of £o for the fine lattice, two for the others). The errors of 
these results are therefore quite small; they were used to investigate the question of when 
£o(£) determined from different S(x,y) reaches its plateau as x and y increase. If we can 
show that the plateau is reached quite early, we can then confidently determine it with 
less statistics in other cases. 

First of all, we should point out the obvious, namely that the existence of a "plateau" 
depends on the accuracy requested; if one demands infinite precision, the plateau would 
of course only be reached as x,y — > oo. Our aim is to determine the relation between £ 
and £o to about 1%. Let us now investigate when a plateau of this accuracy is reached. 

In figures Q — ||| we show the ratios R ss (x, y) and R s t(x, t) and their "ratio of ratios" 
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Figure 2: As in figure [l] for an intermediate coupling, (3 = 6.8, ^ ~ 2. 
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Figure 3: As in figure I for a coarse lattice, (3 = 5.8, £ sa 3. 
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Figure 4: £o estimated from different (ratios of) Wilson loops plotted versus the area 
of the Wilson loops. We show results for a very fine (left) and an intermediate (right) 
lattice. Separate symbols are used to denote results from different y. The solid lines are 
used to indicate our final estimate of £o an d its error. The short-dashed lines denote 
a 1% band around the estimated £o- Points that would otherwise overlap have been 
slightly separated. For details see the main text. 
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Figure 5: As in figure for a coarse lattice, = 5.8, £ = 3. 
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Table 1: Simulation results for the renormalization of the anisotropy, ?7 = £/£o- 
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R ss (x,y)/R s t(x,t) as a function of x for two different values of y. Clearly, excited- 
state contributions cancel in the ratio of ratios long before the individual ratios reach 
their plateaux in x (in fact, it is not clear if they ever do for the unsmeared Wilson 
loops and volumes we use). This is in accord with expectations discussed earlier. More 
quantitatively, it seems that 5(x, y) reaches its fixed-y plateau to high precision when 
x > max(y— 1, 2). With somewhat less accuracy the plateau is apparently already reached 
for x = 2, even for y > 4, at least for not too fine lattices. 

Concerning the plateau of £o(£) consider figures ^| and ||, where the values determined 
from different Wilson loops are presented. For fixed y we see the plateau emerging for 
x > max(y — 1,2). It seems that y = 3 is sufficient to reach the asymptotic plateau of £o 
within the accuracy of interest. In fact, these and other results indicate that the plateau 
can be obtained to 1% or better from just 5(2,3) and <5(3,3).^] In all cases the value so 
obtained is consistent with that from y > 3 and x > y — 1. Of course, at strong coupling 
the errors rise so rapidly as x, y increase, that it is next to impossible to show with 
the same confidence and accuracy as on finer lattices that a plateau has been reached. 
However, within the stated aim of 1% errors our simulations results indicate that this is 
the case. This is also plausible a priori, since x, y = 3, say, which on fine lattices clearly 
gives results within 0.5% or so of the asymptotic value, corresponds to a much larger 
distance in physical units on coarse lattices. 

The final values we quote for £o(£) were obtained by not just considering results from 
fixed x, y (as in the figures above) but also sums of Wilson loops: We sum S(x, y) over 
various x,y satisfying y > 3, x > y — 1 and then find the zero crossing in £q- Larger 
sums were considered for fine lattices, smaller ones for coarse lattices. This was mainly 
done to minimize the subjective element in choosing the final £o(£)- For the errors, on 
the other hand, we were rather conservative and take into account the scatter observed 
between different (sums of) Wilson loops. 

Errors of the individual 5(x,y) or sums thereof (and of £ = t/y in the second 
method) were evaluated with the bootstrap method using 100 bootstrap ensembles. 
Before bootstrapping, the data were binned to check for auto-correlations (and save disk 
space). We always found that errors were stable, at least after some initial binning had 
been performed. Our final results are shown in table [j], where we present the measured 
renormalization of the anisotropy, r) = £/£o- 

Concerning the conceptual underpinnings of our method we should remark that we 
have investigated in one case what happens if one uses smeared Wilson loops. With 
no effort to tune the smearing to be the the same in physical units (however one may 
want to define this) when the (imaginary) heavy quarks are separated along a spatial 
or the temporal direction, we find that the individual ratios R ss (x,y) and R s t(x,t) 
reach their plateaux earlier, as expected, but their ratio reaches its plateau later than 
without smearing. This demonstrates that there is indeed a cancellation of excited-state 
contributions in the "ratio of ratios" of unsmeared Wilson loops. 

5 For fine lattices the £o determined from S(x, 3) for large x > 3 might be just on the boundary of the 
1% band around the asymptotic value we are aiming for. However, it would appear that the £o from 
x = 2, 3 are well within the 1% band, as figures and and all our other results indicate. 
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4 Parameterizing the Renormalization of the Anisotropy 



For future use in simulations we would like to present analytic parameterizations of the 
results obtained in the previous section. Given a good parameterization it is guaranteed 
that observables calculated with the anisotropic Wilson action extrapolate smoothly, 
like a 2 asymptotically, to the continuum limit. If we had performed simulations at the 
central values of the measured £o(£j 9 2 ) of the previous section, data obtained from runs 
with sufficiently high statistics would not lie on a smooth curve. Furthermore, a good 
parameterization allows us to perform simulations at any coupling, not just the ones we 
happened to use in sect. ||. 



4.1 The g 2 dependence of rj 

Obviously, a parameterization should be consistent with known perturbative results. The 
one-loop result for the renormalization of the anisotropy rj of the Wilson gauge action 
was obtained long ago by Karsch |l3| (cf. also [p2|1 ). At one loop rj can be written as, for 
SU(iV), 

V = V(C0,9 2 ) = 1 + ^|^ 2 + 0( 5 4 ) , (4.1) 

with 771 (1) = and 771 (£0 ) increasing monotonically to some finite value in the Hamiltonian 
limit £0 00 (explicit numerical values for the 771 (£0) will be given below). Note that 
to one-loop order we can replace £0 by £ in this equation. We can therefore incorporate 
the known one-loop coefficients 771 into a fit of the results of sect. [|| for fixed £. 

We find that excellent fits of r](g 2 ) for fixed £ are possible to Pade ansatze of the 
form 

2 1 + cig 2 + c 2 g 4 , , 

via = — 7— — 2 — ' ( 4 - 2 ) 
1 + c g 2 

with ci — Co constrained to be consistent with one-loop perturbation theory, i.e. equal 
to 771 (£)/6- The coefficients of parameterizations obtained from fits to such ansatze are 
given in table and a graphic presentation of the £ = 2 and 4 results is exhibited in 
figure |6|. In these and later plots "boosted 1-loop perturbation theory" refers to the 
replacement of the bare coupling g 2 by the "boosted" coupling [^3], El] 

-2 _ y 1 a o\ 

9 ' y/W 8a {l,l)W at {l,l) ' 

A more proper mean-field estimate would require pi] an estimate of the appropriate scale 
for rj, which presumably is much more infra-red than that of the naive boosted coupling 
g 2 . This would lead to better agreement with our non-perturbative determination at 
strong coupling. 



4.2 The £ dependence of 77 

We are also interested in parameterizing our results for rj as a function of £ for fixed j3, 
or, more generally, as a function of both £ and (3. To obtain an idea of how to proceed 
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g 2 =6//3 g 3 =6/|S 



Figure 6: Simulation results and fits of the renormalization of the anisotropy, 77, for fixed 
£ = 2 (left) and £ = 4 (right). For details see the main text. 



£ 


CO 


Cl 


C2 


2.0 


-0.76216 


-0.67686 


-0.04351 


3.0 


-0.79245 


-0.67854 


-0.06625 


4.0 


-0.77715 


-0.64925 


-0.07038 



Table 2: Coefficients of the parameterizations of 77 at fixed £, eq. (|4.2j). 



let us consider what perturbation theory has to say about the £ dependence. A plot of 
£ versus £0 at one-loop using (JO]) and the numerical values of 771 (£0) from [ 13 1 and |22| 
reveals that this relation is amazingly linear, even at large coupling. Rewriting (|4.l| ) as 



£ = ^0 = £o + (£o-i)^5 2 + 0( 5 4 ), fi^^-^-miHo), (4.4) 

shows that this is equivalent to a very weak £0 dependence of fji (£0). This is indeed the 
case, as can be seen in table ^. 

The main £0 dependence of the one-loop result is taken care of by rewriting 771 (£0) in 
terms of 771 (£0). We have tried to parameterize the remainder by fitting 7)1 (£0) to a Pade 
ansatz. We find that the results in table ^ can be represented by 

1.002503 + 0.39100 7/ + 1.47130y 2 - 0.19231 y 3 _ j_ 

V ' " 1 + 0.26287 77 + 1.59008 y 2 - 0.18224 y 3 ' V = £0" ' ' " J) 

This curve reproduces the numbers in table ||| with an accuracy of about 10 -6 , which 
should be sufficient for all practical purposes. 
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su 


'/I \SU / 


1.00 


1.00068(6) 


1.25 


1.01116256 


1.50 


1.01758914 


1.75 


1.02143638 


2.00 


1.02364467 


2.25 


1.02480746 


3.00 


1.02516137 


4.00 


1.02319237 


5.00 


1.02090024 


6.00 


1.01886792 


7.00 


1.01716124 


8.00 


1.01574032 


10.00 


1.01354821 


20.00 


1.00846625 


oo 


1.00250290 



Table 3: One-loop coefficient of rj for the Wilson gauge action (see eqs. (|4.4j ) and (|4,1| ) for 
the precise definition of fji). The £o = 1 result is from |l3[| . All other values are from |22|| ; 
they should be accurate to the precision shown. 



We can now try to find a representation of all data for rj. The most naive hope would 
be that an ansatz of the form 

V U 6 1 + a g A 

might be adequate. We find that this ansatz is not only adequate but provides a nearly 
perfect representation of all data in table |] for oo = —0.77810 and aq = —0.55055 (the 
confidence level Q of a fit is above 0.99, suggesting that our error estimates in sect. ^ 
were indeed rather conservative). It is quite remarkable that all our simulation results 
can be summarized with just two fit parameters. 

Cross sections of the function (4J3) for fixed £ and fixed f3 can be found in figures (?) 



and [8|, respectively. Since there seems to be no systematic trend in the tiny deviations 
of the curves from the data, we recommend the use of ( |4.6[) for future applications of 
the anisotropic Wilson action with 1 < £ < 6 (instead of the fixed-£ representations 
discussed in the previous subsection). 

There is one interesting point concerning the £ dependence of eq. (|4.6| ) that we should 
mention. Recall that in one-loop perturbation theory £ is a linear function of £o to high 
accuracy. The same is therefore true for £o as a function of £. One might think that 
the accuracy of the ansatz ( |4.6j ) means that the same is true beyond the one-loop level. 
However, this is not the case. Trying to fit £o as a linear function of £ for fixed j3 usually 
leads to very small confidence levels (as low as 0.0005) for (5 < 8.0, whereas quadratic fits 
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always work very well.[] The non-linearity can be seen in figure ||| — at least with the help 
of a ruler — where we plot £n versus £ for several (3. This might sound paradoxical at 
first, but the reason is that £ appears on the right hand side of (|4.6|) , not £o; which might 
be considered the natural variable in perturbation theory (though it does not matter on 
the one- loop level). 

First of all, the £ or £o dependence of rji is not the problem; it is much too weak. 
The ansatz ([O]) works just as well if we replace fji by 1.02, say. So we can write ( |4.6| ) as 
77(£) = 1 + bz , where z= (£ — l)/£ to a good approximation, and b = b(g). This implies 
£o = £,/ri = £{1 — bz + b 2 z 2 + . . .). Our above finding concerning fits of £o(£) therefore 
essentially means this: a "quadratic term" (cf. the previous footnote) is present, but 
within errors it is indistinguishable from the term b 2 z 2 that arises from expanding l/rj. 



It is therefore simpler and much better to fit the £ dependence of r\ as in (4J3) instead of 
fitting £o(£) to a quadratic ansatz (equivalently, one can fit £o(£) to the ansatz £o = C/ 7 ? 
with r\ as above). 

To summarize, the non-perturbative £ dependence of rj is within our errors simply 
given by replacing £o by £ in the one-loop formula. Why this should be, we do 
not know, but we are certainly happy to exploit this fact in presenting the simple 
parameterization Q4.6|) of all our data. 



4.3 The derivative <9£o/<9£ 

Among the immediate applications of this result we would like to point out just one, 
related to thermodynamic studies. Obtaining all thermodynamic information requires 
one to take independent derivatives of the partition function with respect to temperature 
and volume. Even if one is ultimately interested the isotropic case, the most natural way 
to take these derivatives is to introduce an anisotropic lattice with independent temporal 
and spatial lattice spacings at an intermediate stage. In addition to the /3-function for 
the dependence of the lattice spacing on the coupling, one needs d^o/d^ to calculate all 
thermodynamic quantities (for details on lattice thermodynamics see e.g. [pq, |26|). 



From eq. (4.6) one obtains 



§({ = !,/) = ! - Mai±i^ = 1 — 7(3) 3 2 . (47) 
at, 6 1 + aog z 

Our 7(5) is shown in figure [l(] together with results from a different non-perturbative 
determination using the "integration of plaquettes" technique (we took the data 

from table 3 in p6| ). This technique can not be used on very coarse lattices and becomes 
very costly on fine lattices. We therefore only have results for 5.7 < (3 < 7.2 to compare 
with. No errors are quoted in for these results, but we tried to estimate them from 
figure 1 in [jl7j (the same information can be found with higher resolution in figure 5.3 
of ref. for example). This error estimate is included in figure [l^ for a few points. 
For our results we also include errors, which were obtained conservatively from fits of r\ 

6 To insure that the ansatz for r/ is well-behaved in the Hamiltonian limit, one might prefer to mean 
£0 = £ + — 1) + 62 (£ — 1) 2 /C by "quadratic ansatz" , instead of the naive £0 = £ + — 1) +&2(£ — l) 2 . 
However, this is not important here; either ansatz works fine for our limited range of £. 
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2 = 6//S 



Figure 7: Fixed-£ cross-sections of our simulation results and the global fit ( f4.6|) for 
renormalization of the anisotropy, ?7 = £/£o- We show our fit up to /? = 5.4. 
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Figure 8: As in figure |7[ but now for fixed-/? cross sections. For fixed (3 the 
renormalization of the anisotropy rj appears to be an essentially linear function of 1 — 1/£, 
even on coarse lattices. 
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Figure 9: Similar to figure g but now we plot the bare versus the renormalized anisotropy 
(for clarity we only show a subset of the (3 values we simulated) . In this case deviations 
from linearity are significant within our errors. For details see the main text. 
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Figure 10: Comparison of non-perturbative determinations of 7 = (1 — d^o/d^)/g 2 at 
£ = 1. The solid line denotes our result, eq. (4/7), for > 5.4, the short-dashed line is 
the result from the integration of plaquette technique. The errors of the two methods 
are shown at several points. We also show results (x) from a previous application of (a 
version of the) ratio method, which gave rather different results. 
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versus £ at fixed (5, instead of from our global fit (|4.6| ) (to be sure, the central values of 
our points in figure are from fl4.7| )). 

Considering that different non-perturbative methods can differ by 0(a 2 ) lattice 
artifacts, the agreement between the two methods in figure 1C is excellent. On the other 
hand, neither of these results agrees very well with estimates from a different application 
of the ratio method [[17], [l8|] , which is also shown in figure For the largest j3 considered 
these results also appear to be in conflict with boosted perturbation theory. (Our results 
here and other results B [3|, ||j have so far always shown boosted perturbation theory to 
agree quite well with non-perturbative determinations for /3 > 7 or so.) The reason for 
this disagreement is not completely clear .[] 

We should mention that in the SU(2) case the integration of plaquette technique (27]] 
and a different non-perturbative method |2^] also agree reasonably well (the latter giving 
slightly larger 7(5)), and are far above the perturbative result. Qualitatively the situation 
is therefore the same as with our SU(3) results. 

Using (^7j) and the /3-function that can be obtained from the accurate scale 
determinations for the isotropic Wilson gauge action in |^9| (and references therein), 
all ingredients required for high precision thermodynamic studies with this action are 
now known for all lattice spacings. Once the string tension and/or Sommer scale ro 
has been determined for the anisotropic Wilson actions, the same will hold for these 
actions. Lattice artifacts in the thermodynamics of the anisotropic Wilson (gauge or 
quark) action are considerably smaller 31 1 than for the isotropic case, so it might be 
interesting to pursue such studies on anisotropic lattices. 



5 Conclusion and Outlook 

Using a method based on ratios of Wilson loops we have presented an accurate, non- 
perturbative determination of the relation between the bare (£0) arid renormalized (£) 
anisotropies of the Wilson gauge action. We have argued that these ratios should be 
thought of as the "finite- volume static potential with excited-state contributions". By 
creating the same physical situation with the heavy quarks separated either along a 
spatial or a temporal direction there should be no significant finite-volume or excited- 
state corrections to the relation between £0 an d £ even for moderately small Wilson loops. 
We have explicitly seen this in our simulations. This fact makes our method quite cheap. 
Nothing in our method is specific to the Wilson action, so it can also be be applied to 
improved gauge actions. 

7 The most likely reason is this: In Jl?] , |l8| the ratios R ss (x, y) and R 3 t(x, £y) were not exactly matched, 
but rather were allowed to differ by an overall factor (the same for all ratios) that was then fitted and 
typically came out a few percent above 1 (unfortunately no errors were quoted for this factor in |l8|). 
From our data we can extract the effect that the deviation of this factor (our "ratio of ratios") from 
its correct asymptotic value 1 has on the final r](£,g 2 ). It turns out that it is sufficient to explain the 
quantitative differences between our and the results in jl?], In particular, for fine lattices f?(£,<7 2 ) 
becomes rapidly more sensitive to deviations of the ratio of ratios from 1, providing an explanation for 
the fact that the difference between our and the results of jl?], Q in figure [To] increases rather than 
decreases for small coupling, contrary to naive expectations. 
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Our results for £,o(£,,g 2 ) are significantly more accurate than previous ones, were 
obtained at a fraction of the cost, and are the only ones that cover the full range 
from weak to strong coupling^] and a large range of anisotropies. Furthermore, we have 
presented a simple parameterization of our results, eq. ( |4.6| ) with ao = —0.77810 and 
a\ = —0.55055, that reproduces all data within errors and is consistent with perturbation 
theory at weak coupling. For lattice spacings of interest in practice we find that the 
renormalization of the anisotropy is much larger than predicted by one-loop perturbation 
theory (even when "boosted"). Our results are therefore crucial for future applications. 

The errors of our results, about 1% on coarse and less on finer lattices, should also 
be sufficient for future applications. Since our parameterization incorporates the known 
one-loop behavior, it is quite clear that the error in an observable due to the remaining 
small error in £o(£> 9 2 ) w ih (almost completely) extrapolate away in the continuum limit|] 

Given these results, the anisotropic Wilson gauge action is as simple to use as the 
isotropic one. The only thing that is currently missing are accurate determinations 
of the physical scale (string tension and Sommer scale ro) of the anisotropic actions. 
This gap should soon be filled, since once the renormalized anisotropy is known, the 
"regular" static potential (cf. sect. [2j) can be determined very accurately. Work on this 
is in progress. 

A number of projects using the anisotropic Wilson action can now start immediately. 
One example is the study of glueballs in pure gauge theory. It will be interesting to 
see how much of the improvement seen in glueball studies with anisotropic improved 
actions |l^] is due to the anisotropy and how much due to the elimination of (most 
of) the O(dg) errors. Another example, as remarked at the end of sect. 4, would be to 
reconsider high precision thermodynamic studies using our results. 

One could also start simulations of heavy quark systems in the quenched 
approximation. An immediate improvement compared to isotropic studies should be 
apparent (cf. p0| ). However, as mentioned in the introduction, truly reliable and 
accurate studies of such systems will have to await the non-perturbative 0(a°) and 
0(a) improvement of Wilson- type quark actions on anisotropic lattices: The first (and 
easier) step is to tune the bare velocity of light of the quarks so that the fermion and 
gauge sectors agree on the renormalized anisotropy. The second step involves the tuning 
of the temporal and spatial clover coefficients to eliminate O(a) errors fllf-FI 

Finally, since our method is free of most lattice artifacts whose elimination would 
be very expensive beyond the quenched approximation, the prospects for extending this 
work to full QCD look good. Of course, simulations for full QCD will be significantly 
more expensive, not just because the determinant of the quark matrix has to be 

Preliminary results indicate that the lattice spacing for (3 = 5.5 and £ > 2 is close to 0.3 fm. The 
range of couplings we covered should therefore be sufficient for all future applications. 

9 We are here referring to the statistical and systematic errors of our method of determining £o(£,(? 2 ); 
another legitimate method can of course differ by 0(a 2 ) terms, which are a priori known to extrapolate 
away in the continuum limit. 

10 This presumably has some (small) effect on the bare velocity of light, which therefore has to be 
retuned iteratively with the clover coefficients. We expect this iterative retuning to converge rapidly, if 
necessary at all. 
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calculated, but also because one has to tune more bare parameters simultaneously to 
obtain consistent quark and gauge anisotropics (now the quark parameters feed back 
into the gauge sector). However, the required tuning might not be as hard as it first 
sounds. 
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